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Abstract 

In the partition function of the Kondo lattice (KL), spin matri- 
ces are exactly replaced by bilinear combinations of Fermi operators 
with the purely imaginary chemical potential A = —mT/2 (Popov 
representation). This new representation of spin operators allows one 
to introduce new Green Functions (GF) with Matsubara frequencies 
u> n = 2irT(n + 1/4) for S = 1/2. A simple temperature diagram tech- 
nique is constructed with the path integral method. This technique 
is standard and does not contain the complicated combinatoric rules 
characteristic of most of the known variants of the diagram techniques 
for spin systems. The effective action for the almost antiferromagnetic 
KL problem is derived. 
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Many systems in statistical physics are described by Hamiltonians con- 
taining spin matrices. Unfortunately, the diagrammatic perturbation the- 
ory for spin systems is complicated. Many variants are based on different 
representations of the spin matrices by Bose or Fermi operators. However, 
unphysical states always arise leading to constraints and complication of the 
Feynman codex. In this paper we construct a simple diagrammatic technique 
(DT) for spin- 1 that differs from the known techniques in the form of the 
GF, but which is standard in other respects, does not contain the compli- 
cated combinatoric rules characteristic of spin systems and permits one to 
take into account the constraints rigorously. 

It is indeed possible to replace exactly spin-| matrices by bilinear combi- 
nations of Fermi operators: 

a z fi -> S z fi = ^(a^fliT - aja u ) ; cr+ -> S+ = a+a u ; aj { -> S n = a+a iT , (1) 

by the basic formula shown in Ref. |TJ (see also Ref. 0]): 

Z = Spe-* ff = ^Spe~^ + i 7V) , (2) 

where Hf is obtained from H by the replacement ([I]), and N = a^ n a\ a . 

icr 

There is no constraint but the purely imaginary chemical potential of pseud- 
ofermions A = —mT/2 leads to the mutual cancellation of the unphysical 
states. 

We analyze here the KL model which is a periodic lattice of magnetic 
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atoms modeled by /-orbitals in a metallic background: 
Ukl = -T t ( t »+^)K^ + JsfT l K^S fl + gY:(H + he^)S' fl . (3) 



We add a uniform (H) and a staggered (h) magnetic field (g = \ib9li where 
fiB is the Bohr magneton and ql is the Lande factor). We consider a simple 
cubic lattice with the notation Q = Qaf = (tt, 7r, tt). Using Popov represen- 
tation of spins, the ratio of the partition function of the interacting system 
to the partition function of the corresponding free system can be represented 
in the form of functional integrals as follows: 
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where the Euclidean action for the KL model is 

S = I dr{J2[^i,ad T ^i,a + a i;a d T a i>a \ - H KL (r)}. (5) 
Jo Y 

We note by Dv the integration over the anticommuting Grassmann variables 
"$? a ,a a . By making the replacement a- ha {r) — * a f( T ) ex P(f| T ) e ^ c > which 
cancels the last term in both exponents in numerator and denominator of 
Eq. (Q), we come to the following boundary conditions for Grassmann fields: 
* CT (/3) = -^(0), = -f CT (0), a a (j3) = -ia a (0), a a {(3) = ia a (0). 

Going over to the momentum representation for all Grassmann variables and 

assuming s sk = ^f p+k a^/ p we obtain: 

p 

s = E^AVEW^^E^S/^+^E^. ( 6 ) 



where the inverse GF of \I/-fields is G 1 = i2iiT(n + 1/2) — + \i with 
dispersion = — E ti,i+se lSk and the inverse GF of a a Grassmann fields is 
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Q^ a = i2nT(m+l/4:) — l/2gHa^ a with unusual Matsubara frequencies. Note 
that Popov representation can be used for spins S — 1 as well. In this case 
the frequencies are shifted to uj m = 27rT(m + 1/3). Moreover the method has 
been also extended to arbitrary spin in Ref. ||. 

We now confine ourselves to the limiting case TV ~ TV ~ T [||] assuming 
the same energy scale for antiferromagnetic (AF) and Kondo correlations. It 
allows us to integrate over the fast \l/ fields with energies e ~ n ^> T using 
the bare electrons GF. We can also integrate over the fast fields a a (u ^> T ) 
taking into account one-site Kondo renormalization of vertices (J s f — > J s f) 
and self-energy parts (G — > G) As a result, a simple DT is constructed. 
Contrary to other DT (see, e.g. 0-0) the constraint on the spin subsystem 
is taken into account exactly. The new action which is written in terms of 
slow ^ and a a variables contains an additional AF Heisenberg interaction 
between spins due to the indirect RKKY exchange 0-0: 

Seff = e g- 1 ^ + e j sf sfrsf - k + s H m 

The last term in (^) can be analyzed separately and represented by auxiliary 
three-component Bose fields 7 (/c) 0: 

SE= E ^Q-^+tX^+Q^ 

fclfc2fT k 

(8) 
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with the following notation: S = -1/4 E( r fc™TVfc0l fc , V* k = (Vk)* and 

k 

Vk = (4>k — /2. In the case Tk <C 2jv only magnetic terms in the effective 
action are important. We note by W the matrix of the quadratic form in 
a a variables. Integrating over all a a fields one can find the nonpolynomial 
action of the AF Heisenberg model in terms of Bose fields 0: 

S H = S [<j) 3 ,v} +logdet^[0 3 ,7?]. (9) 

In the case considered, namely T/v ~ TW ~ T the procedure of derivation of 
the effective action is a little bit complicated. Taking into account the second 
term in (|7|) one has to replace 7 — > 7 — 2j7 s /s^ /ow in (|). As a result, the 
effective action can be rewritten as follows: 

Seff = E n l 7 G ^t7 + 5 o[0 3 , V) + log det W[<f k - 2J sf sl s k l0W }. (10) 

k 

Equations ([7HTD]) are the key result of the present work. This effective action 
describes the low energy properties of the KL model. The last term in (^) 
takes into account the mutual influence of conduction electrons and spins. 
Magnetic instabilities of both kind of electrons could then be easely analysed. 

Let us concentrate on the former problem @. The spin subsystem un- 
dergoes a phase transition with T c = T corresponding to the appearance of a 
nonzero staggered magnetization p as h — > 0. This problem is related to the 
Bose-condensation of the field 0| = + p(/3iV)£5k,Q<L and in one-loop ap- 
proximation results in the usual mean-field equation for AF order parameter 
[0] in the presence of Kondo-scattering processes ||]. Note that a magnetic 
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transition in the localized system may induce a magnetic transition in the 
itinerant system. 

Taking into account the compensation equation 0,0] and calculating 
the logdet iy[0 3 , rj] approximately by the method of stationary phase the 
following expression for the spin subsystem effective action can be obtained: 

s e J f = EfcV + E^xrV 3 -, - i/4£(r™y W-*, (ii) 

where \t and \i are transverse and longitudinal susceptibilities respectively. 
As usual, the transverse susceptibility describes the AF magnons excitations. 
At the temperature range T > T when the condensate solution is absent the 
effective action has the same form except that the transverse and longitudinal 
susceptibilities describe the paramagnon excitations which can result in some 
untrivial effects in heavy- fermion compounds ||. These excitations introduce 
a new energy scale corresponding to the critical behaviour. 

Summarizing, we constructed a simple diagrammatic technique which 
allows one to analyze the effective action of the KL model when the energy 
scales for AF and Kondo correlations are the same. This effective action 
describes the slow electron subsystem interacting with the spin fluctuations 
of either magnon or paramagnon type. 
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